Modeling Human Dynamics with Adaptive Interest 
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Recently, increasing empirical evidence indicates the extensive existence of heavy tails in the 
interevent time distributions of various human behaviors. Based on the queuing theory, the Barabasi 
model and its variations suggest the highest-priority-first protocol a potential origin of those heavy 
tails. However, some human activity patterns, also displaying the heavy-tailed temporal statistics, 
could not be explained by a task-based mechanism. In this paper, different from the mainstream, we 
propose an interest-based model. Both the simulation and analysis indicate a power-law interevent 
time distribution with exponent -1, which is in accordance with some empirical observations in 
human- initiated systems. 
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I. INTRODUCTION 

Human behavior, as an academic issue in science, has a 
history of about one century from Watson Q. As a joint 
interest of sociology, psychology and economics, human 
behavior has been extensively investigated during the last 
decades. However, due to the complexity and diversity 
of our behaviors, the in-depth understanding of human 
activities is still a long-standing challenge thus far. Actu- 
ally, in most of the previous works, the individual activ- 
ity pattern is usually simplified as a completely random 
point-process, which can be well described by the Poisson 
process, leading to an exponential interevent time distri- 
bution That is to say, the time difference between 
two consecutive events should be almost uniform, and 
the long gap is hardly to be observed. However, recent! 
the empirical studies on e-mail 0] and surface mail 
communication show a far different scenario: those com- 
munication patterns follow non-Poisson statistics, char- 
acterized by bursts of rapidly occurring events separated 
by long gaps. Correspondingly, the interevent time dis- 
tribution has a much heavier tail than the one predicted 
by an exponential distribution. The heavy tails have also 
been observed in many other human behaviors S @| i in- 
cluding market transaction 0, Q , web browsing , movie 
watching short message sending and so on. The 
increasing evidence of non-Poisson statistics of human ac- 
tivity pattern highlights a question: what is the origin of 
those heavy tails? Based on the queuing theory, Barabasi 
et al. proposed a simple model [1, [H, where the indi- 
vidual executes the highest-priority task first, and they 
suggested the highest-priority-first (HPF) protocol a po- 
tential origin of those heavy tails. 

The queuing model gets a great success in explaining 
the heavy tails in many human-oriented dynamics. How- 
ever, some other human activity patterns, also display- 
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ing the similar heavy-tailed phenomenon, could not be 
explained by a task-based mechanism. For example, the 
actions on browsing webs watching on-line movies 
[lo| . and playing on-line games are mainly driven 
by the personal interests, which could not be treated 
as tasks needing to be executed. The in-depth under- 
standing of the non-Poisson statistics in those interest- 
driven systems requires a new model out of the perspec- 
tive of queuing theory. In this paper, different from the 
mainstream task-based models, we propose an interest- 
based model. Both the simulation and analysis indicate 
a power-law interevent time distribution with exponent 
-1, which is in accordance with some empirical human- 
initiated systems. 



II. MODEL 

Before introducing the mathematical rules of our 
model, let us think of the changing process of our in- 
terests on web browsing according to our daily experi- 
ences. If a person has a long period not browsing the 
web, an accidental visit may give him a good feeling and 
wake his interest on the web browsing. Next, during the 
actions, the good feeling continues and the frequency of 
web browsing may increase. Then, if the frequency is too 
high, he may worry about it, thus reduces those browsing 
actions. Such similar experiences can be found in many 
other daily actions, such as playing games, seeing movies, 
and so on. In a word, we usually adjust the frequency of 
the daily actions according to our interest: greater inter- 
est will lead to higher frequency, and vice versa. Some 
simple assumptions extracted from our daily experiences 
are as follows: Firstly, for a given interest-driven behav- 
ior, each action will change the current interest, while 
the frequency of actions depends on the interest. It likes 
an active walker [l^, [l^ , whose motion is affected by the 
energy landscape, while the motion track could simulta- 
neously change the landscape. Secondly, the interevent 
time T has two thresholds: when r is too small (i.e.. 
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FIG. 1: (upper panel) The succession of events predicted by 
the present model. The total number of events shown here 
is 375 during 10® time steps, (lower panel) The correspond- 
ing changes of r{t). The data points are obtained with the 
parameters ao = 0.5 and T2 — 10*. 

events happen too frequently), the interest will be de- 
pressed, thus the interevent time will increase; while if 
the time gap is too long, we will increase the interest to 
mimic its resuscitation induced by a casual action. 

According to those assumptions, we propose an 
interest-based model as follows: (i) The time is discrete 
and labelled by t = 0, 1, 2, • ■ • , the occurring probabil- 
ity of an event at time step t is denoted by r{t). The 
time interval between two consecutive events is call the 
interevent time and denoted by r. (ii) If the {i + l)th 
event occurred at time step i, the value of r is updated 
as r{t + 1) = a{t)r{t), where 

{ao, n < Ti, 

ao\ >T2, (1) 

a(t-l), Ti <n <T2. 

If no event occurred at time step t, we set a(t) — a{t — 
1), namely a{t) keeps unchanged. In this definition, Ti 
and T2 are two thresholds satisfied Ti <C T2, Ti is the 
time interval between the (i -I- l)th and the ith events, 
and ao is a parameter controlling the changing rate of 
occurrence probability (0 < ao < 1). If no event happens, 
r will not change. Clearly, simultaneously enlarge (by the 
same multiple) Ti, T2, and the minimal perceptible time, 
the statistics of this system will not change. Therefore, 
without lose of generality, we set Ti = 1. 



III. SIMULATION AND ANALYSIS 

In the simulations, the initial value of r is set as 
''0 — f{i = 0) = which is also the possibly max- 
imal value of r{t) in the whole simulation process. As 
shown in Fig. 1, the succession of events predicted by 
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FIG. 2: (Color online) The interevent time distributions in 
log- log plots, (a) Given ao = 0.5, Pij) for different T2, where 
the black, red and green curves denote the cases of T2 = 
10^ 10^ and lO", respectively, (b) Given T2 = lO", P(r) for 
different ao, where the black, red and green curves denote the 
cases of ao = 0.8, 0.5 and 0.2, respectively. The black dash 
lines in both (a) and (b) have slope -1. Each distribution 
contains 10® events. 

the present model exhibits very long inactive periods 
that separate the bursts of rapidly occurring events, and 
the corresponding r{t) shows a clearly seasonal property 
(quasi-periodic behavior). Actually, in a period, the max- 
imal and minimal values of r[t) respectively determined 
by Ti and T2 as r^ax ~ and r^i_-a ^ ?2^^- This quasi- 
periodic property will be applied in the further analysis. 
Note that, in a specific quasi-period, r^ax can be smaller 
than T-f ^ and rmin can be smaller than T^^. It is because 
T < Ti could happen when r{t) < T^^ and r < T2 could 
happen when r(t) < T2~^. 

Figure 2 reports the simulation results with tunable T2 
and ao. Given ao = 0.5, if T2 ^ Ti, the interevent time 
distribution generated by the present model displays a 
clearly power law with exponent -1; while if T2 is not suf- 
ficiently large, the distribution P{t) exhibits a departure 
from a power-law form with a cut-off in its tail. Corre- 
spondingly, given sufhciently large T2, the effect of ao is 
very slight, thus can be ignored. 

Taking into account the quasi-periodic property of r(t), 
we raise two approximated assumptions before analyti- 
cal derivation: (i) The statistical property of P{t) is the 
same as that in a single period; (ii) Within one period, 
the statistical property of P{t) in the r-increasing half is 
the same as that in the r-decreasing half. In the reducing 
process, r{t) = rmOQ, where z = 0, 1, 2, • ■ •, I. The integer 
/ denotes the number of the events in the reducing pro- 
cess (also the number of different values of r{t)), whose 
value is about 

/«-log„jT2/Ti) (2) 

since r^ax ^ T^^^ and rmin ^ T2^- fm is the initial value 
(it is also the maximum value) of r{t) in a reducing pro- 
cess. Note that, for different reducing processes, the val- 
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ues of rm are not always the same. Though has the 
same order of magnitude with T^^ = 1.0, its value can 
be less than T^^ in a specific process. The average value 
of will be calculated later in this paper. 

If the current occurring probability is r{t) — rmO-o, the 
probability that the next event will happen at the time 
t + Tis: 

Q{t) = (1 - r.^al,y~\^al (3) 

Considering every value of r{t) in the reducing process, 
the interevent time distribution of the reducing process 
is: 

I 

P(r)=/-i^(l-r™a^r-V™4. (4) 

i=0 

According to the approximated assumptions above, the 
interevent time distribution of all the successions can also 
be expressed by Eq. (4), which can be approximately 
rewritten in a continuous form, as: 

P(t) « I {I- rraaiy-^rr.-.a'^^dx. (5) 
Jo 

Therefore, P{t) can be further expressed as: 

P(r) « -[(1 - rm4Y - (1 - r.raY]{\na^)-^r\-\ (6) 

From Eq. (6), for a fixed r„i, when / is large enough (it is 
equivalent to the condition T2 3> Ti), P{t) has a power- 
law tail with exponent -1. In addition, this analytical 
result also provides an explanation about the departure 
from a power law when T2 is not sufficiently large. 

As discussed before, for different reducing processes 
of r{t)^ the possible values of are not always the 
same (see also the lower panel of Fig. 1: for different 
quasi-periods, the maximum values of r{t) are different). 
Since the order of magnitude of is comparable with 
Tj~^ = 1.0 (it is equal to rg), the minimum value of 
f{t), r^aQ, has the same order of magnitude with roa^. 
Making the approximated assumption that the minimum 
value of r{t) is given by roag in a r-increasing process, 
and the maximum value of r{t) in the next r-decreasing 
process is roa§ (roo^ is also the start point in the next de- 
creasing process) , then the probability density of k reads 

/-fe-i 

^{k)^ro4 n (l-^o«o"'')- (7) 

i=0 

Therefore, the average value of rm is 

/-I /-I /-fc-i 

fc=0 fc=0 1=0 

(8) 

This average value of r™ calculated by Eq. (8), as well as 
the integer part of — log^j^ (T2/T1) (as the approximation 
of /), can be directly used in the approximate calculations 




10° 10' 10- 10' 10' 10' 10' 



FIG. 3: The comparison of the analytical (black solid line) and 
numerical (gray circles) results of interevent time distribution. 
The numerical data are obtained with parameters ro = 1.0, 
ao = 0.5 and T2 = 10*. The analytical results are calculated 
by Eq. (6) with ao = 0.5, 7 = 13 and r„ = 0.50. The 
black dash line has slope -1. The numerical results contain 
lO'' events. 



of Eq. (6). Given ro = 1.0, ao = 0.5, T2 = 10^ and Ti = 
1, one obtains / « - log„^(T2/Ti) = 13, and (r^) « 0.50 
by Eq. (8). Accordingly, Fig. 3 reports the comparison 
of analytical and simulation results, which are well in 
accordance with each other. 



IV. CONCLUSION AND DISCUSSION 

A novel model on human dynamics is proposed in this 
paper. Different from the mainstream queuing models, 
the current model is driven by the personal interests. In 
this model, the frequency of events are determined by the 
interest, while the interest are simultaneously affected by 
the occurrence of events. This interplay working mech- 
anism, similar to the active walk |15l . Il6| . is a genetic 
origin of complexity of many real-life systems. The rules 
in the current model are extracted from our daily life, 
and both the analytical and simulation results agree well 
with the empirical observation, such as the activity pat- 
tern of web browsing Q. Our work indicates a simple 
and universal mechanism in human dynamics, that is, a 
people could adaptively adjust their interest on a specific 
behavior (e.g. watching TV, browsing web, playing on- 
line game, etc.), which leads to a quasi-periodic change 
of interest, and this quasi-periodic property eventually 
gives raise to the departure of Poisson statistics. 

Besides the HPF protocol and the current model, there 
are also some other mechanisms that can lead to a power- 
law interevent time distribution. For example. Hidalgo 
[l3| pointed out that a Poissonian individual with char- 
acteristic time varying randomly in time could generate 
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a power-law interevent time distribution with exponent 
-2. In addition, Vazquez [l^ showed that if the current 
executing rate is hnearly correlated with the average ex- 
ecuting rate in a immediate predecessor period, the in- 
terevent time distribution will follow a power-law form. 

Note that, although in the recent empirical works, the 
power-law form is widely used to fit the interevent time 
distribution of human behaviors, there exists a debate 
about the choice of fitting functions for this distribution 
in the e-mail communication [Tol . [20| . Actually, a candi- 
date, namely log-normal distribution, has also been sug- 
gested [l^ to describe the non-Poisson temporal statis- 
tics of human activities. The stretched exponential dis- 
tribution [21I, , interpolating between a power law and 
an exponential form, serves as another candidate (see, for 
example, the distribution of interevent time between two 
consecutive transactions initiated by a stock broker p^). 
A clear understanding of the tails in the interevent time 
distribution asks for in-depth exploration on empirical 
data in the future. 

The concept and methodologies related to the statis- 
tics of the interevent time can also find its applications in 
some other systems. For example, similar statistical anal- 
ysis can be addressed on the spacing between the consec- 
utive occurrences of the same letter in written text 
and the time difference between successive events above 
a certain threshold (i.e., extreme events) [23|. 

Finally, we point out some limitations in the current 
model. Firstly, it can only generate the power-law in- 
terevent time distribution with exponent -1, which does 



not agree with some real human-initiated systems with 
different power-law exponents. Secondly, we assume that 
the changing rate of the occurring probability, oq, is fixed 
as a constant in every rising or decaying process. This 
assumption is very ideal, and we could not find any sup- 
port from the empirical data. Third, as stated by Kentsis 
[2J|, there are countless ingredients affecting the human 
dynamics, and for most of them, we do not know their im- 
pacts. Those ingredients, such as the social content, the 
semantic content and the periodicity due to circadian and 
weekly cycles, have not been considered in the present 
model, neither the HPF protocol. However, although 
this model is rough and may contain some artificial as- 
sumptions, it provides a start point of modeling interest- 
based human dynamics. The human-initiated systems 
are the most complex systems, and there must be many 
underlying mechanisms having not been discovered yet. 
We believe our model could highlight the readers in this 
rapidly growing area. 
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